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ABSTRACT 


This  report  presents  results  of  laboratory  calibrations  of  the  EO-1  Advanced  Land  Imager 
(ALI)  instrument.  Specifically,  it  is  concerned  with  the  estimated  spatial  modulation  transfer 
function,  or  MTF.  A  combination  of  design  data  for  the  detectors,  measurements  on  the  optical 
subsystem,  and  laboratory  measurements  of  the  complete  instrument  are  used  to  construct  a  model  of 
the  system  spatial  transfer  function  (STF)  of  the  ALI.  The  optical  subsystem  measurements,  done  by 
the  telescope  builder,  SSG,  consisted  of  interferograms,  which  they  analyzed  with  Zemike 
polynomials.  Those  polynomials  were  used  to  model  the  optical  transfer  function.  The  system 
measurements  were  scans  of  a  knife-edge,  viewed  through  an  imaging  collimator.  The  parameters  of 
the  STF  model,  in  particular  focus  error  and  carrier  diffusion  in  the  detectors,  were  adjusted  to  obtain 
the  best  agreement  with  the  system  measurements.  The  full,  two-dimensional  STFs  for  the  ten 
spectral  bands  of  the  instrument  are  presented. 
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1  INTRODUCTION 


The  Advanced  Land  Imager  (ALI)  is  the  primary  instrument  on  the  Earth  Observing- 1  (EO-1) 
spacecraft  of  the  NASA  New  Millennium  program.  A  general  description  of  the  ALI  has  been  given 
by  Bicknell,  et  al.^  For  an  overview  of  the  instrument’s  calibration,  see  Lencioni,  et  al.^  The  spatial 
response  of  the  ALI  instrument  was  calibrated  in  the  laboratory  during  the  period  July  to  December, 
1998,  prior  to  its  integration  with  the  spacecraft.  Heam,  Mendenhall,  and  Willard^  described 
techniques  used  for  those  calibrations.  This  document  presents  the  final  results  for  the  calibrations  of 
the  Modulation  Transfer  Function  (MTF).  In  the  following  sections,  the  measurements  made  in  the 
laboratory  are  first  summarized,  then  the  steps  in  processing  and  analyzing  those  measurements  are 
outlined.  Details  of  the  MTF  model  developed  for  the  instrument  are  then  presented.  Finally,  results 
for  the  various  bands  are  shown  graphically.  Numerical  values  of  the  MTFs  are  presented  in 
Appendix  A. 

Spatial  response  is  fully  characterized  by  a  spatial  transfer  function  (STF).  The  full  STF  of  the 
instrument  is  a  two-dimensional  complex  function,  which  depends  on  the  wavelength  band,  and  to  a 
lesser  degree,  on  the  position  within  the  field  of  view.  The  STF  is  the  Fourier  transform  of  the  point- 
spread  function  (PSF).  The  modulation  transfer  function  is  normally  defined  as  the  magnitude  of  the 
STF,  a  real,  non-negative  quantity.  An  additional  phase  transfer  function  (PTF)  is  then  required  to 
describe  the  spatial  response  completely. 

In  the  ALI,  as  in  most  good  imaging  instruments,  the  PSF  is  essentially  symmetrical.  As  a 
result,  the  imaginary  part  of  the  STF  (Im(STF))  is  very  small,  compared  to  the  real  part  (Re(STF)). 
The  conventional  MTF  then  looks  like  the  absolute  value  of  Re(STF).  For  simplicity,  we  prefer  to 
work  with  the  STF,  rather  than  the  MTF  and  a  PTF  which  goes  through  sudden  swings  of  1 80°  where 
Re(STF)  changes  sign.  We  have  habitually  called  the  STF  the  “complex  MTF,”  or  simply  the  MTF 
while  working  on  the  ALI  calibration.  In  the  rest  of  this  document,  the  terms  have  been  used 
interchangeably.  The  reader  should  simply  be  aware  that  we  are  not  following  the  conventional 
terminology. 
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2  KNIFE-EDGE  SCANS  PERFORMED 


The  approach  we  used  to  measure  the  MTFs  of  the  ALI  was  to  project  the  image  of  a  knife- 
edge  into  the  instrument,  and  scan  the  edge  slowly  across  a  selected  set  of  detectors  while  the 
instrument  records  the  image  data  at  its  nominal  rates.  The  knife-edge  was  physically  located  at  the 
focus  of  a  collimator  and  the  scan  was  performed  by  translating  its  support  stage  slowly  in  the  plane 
of  best  focus.  For  each  pixel  in  the  set,  the  detected  signal  vs.  time,  after  normalization,  yields  an 
edge-spread  function  (ESF).  Differentiation  of  the  ESF  gives  the  line-spread  function  (LSF).  Fourier 
transformation  of  the  LSF  results  in  the  STF,  or  complex  MTF,  along  one  axis. 

All  measurements  reported  here  were  made  while  the  instrument  was  in  a  thermal  vacuum 
chamber,  at  its  normal  operating  temperatures.  For  each  of  the  four  sensor  chip  assemblies  (SCAs), 
the  knife-edge  was  scanned  across  a  set  of  pixels  near  the  middle  of  that  SCA,  for  each  band.  Scans 
were  performed  in  both  the  cross-track  and  in-track  directions  of  the  instrument  by  using  two  sides  of 
the  9  mm  open  square  in  a  negative  “USAF  1951”  target  as  our  knife-edges.  (The  rest  of  the  target 
was  covered  by  black  plastic.)  When  switching  from  one  SCA  to  the  next,  the  instrument  was  re¬ 
positioned  on  an  azimuthal  rotation  stage  in  the  vacuum  chamber  in  order  to  scan  each  SCA  with  the 
knife-edge  within  -0.5°  of  the  collimator  axis. 

The  knife-edge  was  moved  at  200  pm/s,  which  translates  to  127.4pm/s  at  the  ALI  focal  plane. 
The  frame  rates  were  226  s'‘  for  the  multispectral  (MS)  bands,  and  678  s"’  for  the  panchromatic 
(Pan)  band.  The  spatial  sampling  intervals  were  0.564  pm  and  0.188  pm  for  the  MS  and  Pan  bands, 
respectively,  or  approximately  70  samples  per  pixel  width  in  all  cases. 

The  recorded  data  files  contained  10,  20,  or  30  seconds  of  data  each,  from  the  entire  focal 
plane.  As  many  as  six  of  the  MS  bands  could  be  scanned  at  once  in  the  cross-track  direction.  No  more 
than  three  bands  could  be  scanned  in  one  pass  in  the  in-track  direction.  Approximately  four  dozen 
scan  data  files  were  ultimately  recorded,  once  the  measurement  system  was  deemed  to  be  working 
correctly. 
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3  DATA  PROCESSING 


Analysis  of  the  data  from  the  knife-edge  scans  proceeded  in  several  steps.  While  the  scans  were 
being  performed  in  the  thermal  vacuum  chamber  in  the  clean  room,  a  portion  of  the  data  was 
analyzed  and  graphically  displayed  on  a  monitor  next  to  the  control  computer.  This  Quick-Look 
analysis  was  a  necessary  check  that  the  equipment  was  set  up  and  working  properly.  The  data  files 
were  stored  on  the  hard  drives  of  a  Silicon  Graphics  R- 10000  workstation  (called  the  Performance 
Assessment  Machine  (PAM)),  and  periodically  backed  up  on  DLT  cartridges. 

3.1  PRE-PROCESSING 

After  all  of  the  knife-edge  scan  files  were  recorded,  preliminary  processing  was  done  on  the 
PAM.  A  program  was  run  to  de-scramble  the  raw  data  (“mandescram,”  written  in  C).  Next,  running 
under  IDL*,  the  “get  knife  scan”  procedure  read  the  full  MS  or  Pan  data  file  and  wrote  a  much  more 
compact  file  (*.kes)t,  containing  only  the  data  from  the  SCA  and  bands  of  that  particular  scan. 

For  convenience,  further  analysis  under  IDL  was  performed  after  copying  the  *.kes  files  to  a 
180  MHz  Power  Macintosh.  (It  was  found  to  be  just  as  fast  as  the  PAM  for  running  IDL  procedures, 
since  its  clock  rate  is  the  same.  It  simply  lacked  the  massive  storage  capability  of  the  PAM.) 

The  next  step  was  to  analyze  the  signal  from  each  detector  to  derive  its  edge-spread  function. 
A  sigmoid  curve  (hyperbolic  tangent)  was  fitted  to  each  pixel’s  ESF,  in  order  to  resample  the  ESFs  to 
co-align  all  of  the  edge-crossing  times.  In  the  process,  automatic  selection  criteria  were  applied  to 
reject  pixels  that  had  data  artifacts  or  incomplete  scans.  An  intermediate  file  (*.kes.fit)  was  written 
to  facilitate  the  further  processing  steps.  All  of  these  steps  were  performed  by  IDL  procedure 
“MTF_process  1  .pro” 

3.2  FULL  DATA  ANALYSIS 

To  transform  the  ESFs  to  MTFs,  another  procedure  (“MTF_process8.pro”)  was  developed  to 
read  the  *.kes  and  *.kes.fit  files,  then  derive  the  LSF  and  MTF  functions  from  them.  For  each  SCA 
and  band  combination,  the  complex  MTF  was  computed  on  a  pixel-by-pixel  basis,  then  the  mean  and 
standard  deviation  was  derived.  These  are  written  to  an  IDL  data  file  (*.mtf)  for  further  reference, 
and  plots  are  generated  to  display  the  results  of  each  step  in  the  analysis. 

Finally,  for  each  band  the  complex  MTFs  for  the  individual  SCAs  are  averaged  and  plotted 
together,  along  with  a  curve  representing  the  MTF  model.  Again,  the  results  are  written  to  an  IDL 
data  file. 


♦  Interactive  Data  Language,  from  Research  Systems,  Inc.,  4990  Pearl  East  Circle,  Boulder,  CO 
80301,  (see  www.rsinc.com). 

t  In  this  context,  *  represents  a  “wildcard,”  substituting  for  the  root  name  of  the  data  file. 
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3.3  TYPICAL  MTF  SCAN  AND  ANALYSIS 

A  typical  page  of  plots  generated  in  the  analysis  of  one  scan  file  is  shown  in  Figures  1  to 
4.  In  all  plots  of  this  sort,  the  independent  variable  is  spatial  position  at  the  focal  plane  (in 
pm)  or  frequency  in  image  space  (in  cycles/mm).  These  are  convertible  to  object  space  by 
dividing  or  multiplying  by  the  ALI  focal  length  (0.946  m)  to  obtain  units  of  microradians  or 
cycles/mrad,  respectively. 

3.4  EQUIPMENT  DEFECTS  AND  ARTIFACTS 

Not  all  of  the  scans  produced  data  as  clean  as  that  seen  in  Figures  1  to  4.  Artifacts  were 
sometimes  seen  in  the  data,  particularly  from  in-track  scans,  wherein  the  knife-edge  moved 
vertically.  They  included  occasional  jitter,  with  amplitude  on  the  order  of  a  micrometer,  and 
apparent  errors  in  translation  speed.  It  was  later  discovered  that  the  vertical  Newport  slide  was  not 
properly  adjusted.  Fortunately,  there  is  enough  redundancy  in  the  set  of  scans  performed  that  those 
that  have  obvious  artifacts  could  simply  be  omitted  from  the  MTF  averages. 

There  are  two  “leaky”  pixels  among  the  11,520  MS  pixels.  Their  signals  are  cross-coupled  to 
all  of  the  other  pixels  in  the  same  row  of  their  particular  SCA.  Fortunately,  they  affect  only  the 
even-  or  odd-numbered  pixels  of  that  band/SCA  combination.  It  was  an  easy  matter  to  exclude  the 
contaminated  pixels  from  the  MTF  analyses. 
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Figure  1.  Edge-spread  functions  of  SCA2,  Band  4  pixels  scanned  in  the  cross-track  direction, 
from  file  no.  16638.  Pixels  not  used  in  further  analysis  are  shown  with  dotted  lines. 
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Figure  2.  Average  line-spread  function  derived  from  the  data  shown  in  Figure  1. 
The  dotted  lines  are  one  standard  deviation  above  and  below  the  mean. 
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Figure  3.  Normalized  transfer  functions  derived  from  the  data  shown  in  Figure  1. 
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Figure  4.  Average  normalized  transfer  function  derived  from  the  data  shown  in  Figure  1.  The  dotted  lines 
are  one  standard  deviation  above  and  below  the  mean.  The  curve  labeled  "Diffraction  +  detector’’  is  computed 
from  the  simplest  model  of  diffraction  by  a  circular  aperture  and  a  sharp-edged  detector  response.  The  specified 
MTF  at  the  Nyquist  frequency  is  indicated  with  dotted  lines. 
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4  MTFMODEL 


Our  goal  is  to  describe  the  MTF  of  the  instrument  with  a  model  which  contains  only  enough 
parameters  to  fit  the  measurement  data  within  the  experimental  uncertainties,  and  which  will  enable 
other  investigators  to  compute  the  MTF  for  their  own  needs.  It  is  important  to  have  a  model  for  the 
MTF,  since  the  knife-edge  scan  measurements  can  only  produce  one-dimensional  cross-sections  of 
the  MTF,  which  is  a  two-dimensional  function.  The  following  subsections  describe  the  various  factors 
which,  when  multiplied  together,  yield  the  system  MTF  of  the  ALI  instrument.  These  are  the  optical 
system  MTF,  the  detector  MTF,  and  the  MTF  of  the  scan  motion  of  the  instrument  as  it  collects 
images  on  orbit.  Uncertainties  in  the  model  are  summarized  at  the  end  of  this  section. 

4.1  OPTICAL  MTF 

The  optical  MTF  estimated  from  our  laboratory  measurements  results  from  the  combined 
optical  transfer  functions  (OTFs)  of  the  ALI  telescope  and  the  collimator.  It  is  modified  by  the 
effects  of  light  scattering  from  the  imperfect  surfaces  of  the  ALI  mirrors. 

4.1.1  ALI  telescope 

The  OTF  of  an  optical  system  can  be  obtained  from  the  autocorrelation  of  the  total  wavefront 
error  (WFE)  over  the  pupil  of  the  system.  SSG,  Inc.,  the  builder  of  the  ALI  optical  system,  has  used  a 
laser  unequal-path  interferometer  (LUPI)  to  characterize  the  WFE  of  the  ALI.  They  obtained 
interferograms  at  a  dozen  points  distributed  over  the  1.26°xl5°  field  of  view  of  the  optical  system. 
The  interferograms  were  referenced  to  a  fixture  representing  the  focal  plane.  The  WFE  derived  from 
each  interferogram  was  fitted  with  a  set  of  37  Zemike  polynomials.  The  resulting  coefficients 
encapsulate  the  WFE  estimate.  The  first  three  Zemike  polynomials,  representing  piston,  tip,  and 
tilt,  are  set  to  zero.  The  fourth  coefficient  represents  a  focus  error. 

We  use  the  Zemike  coefficients  from  SSG  to  obtain  an  interpolated  set  of  coefficients  for  any 
given  point  within  the  MS/Pan  array  of  the  ALI.  Since  the  focal  plane  fixture  used  by  SSG  was  not  at 
the  correct  focus  of  the  system,  the  focus  terms  have  to  be  modified  to  represent  the  position  of  the 
ALI  focal  plane  array  (FPA)  after  final  shimming.  Three  parameters,  consisting  of  the  axial  shift 
(piston),  and  two  tilt  angles  control  this  correction. 

To  ascertain  the  amount  of  focus  error  in  the  integrated  system,  MTF  scans  were  performed 
with  varying  amounts  of  defocus  applied  at  the  collimator,  as  described  in  Reference  1.  The  figure  of 
merit  (FoM)  chosen  to  test  the  defocus  was  the  integral  of  the  MTF  from  zero  to  the  detector 
sampling  frequency  of  the  FPA.  A  plot  of  this  FoM  as  a  function  of  the  axial  displacement  of  the 
knife-edge  from  the  collimator  focus  is  shown  in  Figure  5.  The  plot  also  shows  a  parabola  fitted  to 
the  highest  values,  to  derive  the  focus  error.  These  results  indicate  that  the  ALI  FPA,  at  operating 
temperature,  is  approximately  75±20  pm  away  from  the  optimum  focus  position.  This  does  not 
seriously  degrade  the  optical  resolution,  but  must  be  taken  into  account  when  modeling  the  system 
MTF. 
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Figure  5.  ALI  focus  figure  of  merit  vs.  axial  displacement  of  the  knife-edge  from  the  collimator  focus 
(Z  Offset).  Data  symbols  represent  individual  MTF  scans.  The  curve  is  a  parabola  fitted 
to  the  points  within  500  pm  of  the  peak. 


4.1.2  Collimator 

Wavefront  error  introduced  by  the  collimator  should  be  added  to  the  WFE  of  the  ALI  optics 
before  the  PSF  or  MTF  is  calculated.  The  collimator,  and  the  methods  used  to  set  it  to  best  focus,  are 
described  in  a  paper  by  Willard.^  The  spherical  aberration  of  the  collimator  is  calculated  to  be  0.007 
pm  rms,  when  the  collimator  and  ALI  pupils  are  correctly  aligned.  The  remaining  rms  WFE,  caused 
mainly  by  the  beamsplitter,  is  approximately  0.032  pm.  Though  the  spherical  aberration  is  easily 
evaluated,  we  lack  a  full  description  of  the  collimator  WFE. 

4.1.3  Mirror  scatter 

Microscopic  roughness  of  the  surfaces  of  some  of  the  ALI  mirrors  produces  a  large-angle 
scatter  of  a  fraction  of  the  incident  light.^  The  MTF  of  this  scattered  light  is  a  very  narrow  peak 
around  zero  spatial  frequency.  The  total  optical  MTF  is  the  sum  of  this  narrow  scattering  peak  and 
the  specularly-reflected,  normal  MTF.  The  fractional  division  between  the  scattered  and  specular 
light  varies  inversely  with  the  fourth  power  of  the  wavelength.  Our  MTF  model  includes  a  parameter 
to  allow  for  this  scattering  MTF.  It  is  adjusted  to  match  the  MTF  measurements,  though  since  the 
knife-edge  scans  cover  a  limited  angular  range,  they  are  not  sufficient  to  determine  the  fullest  extent 
of  the  scattering.  Scattering  accounts  for  ~5%  of  the  light  at  the  shortest  ALI  wavelengths  (band  1’, 
433-453  nm). 
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4,2  DETECTOR  MTF 


The  spatial  response  of  the  detectors  was  initially  modeled  as  unity  within  a  defining 
rectangular  mask  (39.6  pm  x  40  pm  for  the  MS  pixels,  and  13.2  pm  square  for  the  Pan  pixels),  and 
zero  outside  of  that  area.  The  corresponding  pixel  MTF  is  a  product  of  two  sine  functions: 

=  sinc{wf^)sinc(hfyj 


where  w  and  h  are  the  width  and  height  of  the  detector  mask,  respectively,  and  (f^  fy)  are  the 
components  of  the  spatial  frequency. 


Measured  values  of  the  ALI  system  MTF  tended  to  decrease  approximately  linearly,  in 
comparison  with  the  system  MTF  computed  using  this  simple  detector  model,  after  adjustment  of  the 
free  optical  parameters.  Good  fits  between  measured  and  modeled  MTFs  are  only  obtained  when  a 
carrier  diffusion  MTF  is  included  in  the  model^: 

^TFdet  —  ^TFpjx  ■  MTFej)  , 


MTFed  =  exp 


// 


% 


f  = 


+  /v 


where  f  and  g  are  constants. 

The  VNIR  detectors  (bands  1  ’  through  4’  and  Pan)  are  Si  PIN  photodiodes,  formed  within  the 
readout  circuit  (ROIC).  Good  MTF  fits  were  obtained  with  g  =  \,fo  =  200±50  cycles/mm,  and  the 
nominal  geometric  size  of  the  mask  regions.  Thus  the  spatial  response  of  the  VNIR  detectors  is  a 
convolution  of  the  nominal  mask  rectangle  with  a  Lorentzian  function  having  a  full  width  at  half- 
maximum  (FWHM)  of  approximately  0.8  pm. 

The  HgCdTe  SWIR  detectors  (bands  5’,  5,  and  7)  are  on  a  chip  bump-bonded  to  the  ROIC.  Best 
MTF  fits  were  obtained  in  this  case  with  g  =  1.5,yj)  =  35  cycles/mm,  and  the  geometric  size  reduced  to 
36.8  pm  square. 

A  plot  of  the  two-dimensional  MTF  of  the  Pan  detectors  is  shown  in  Figure  6.  Figure  7  shows 
the  Pan  detector  spatial  response  (point-spread  function).  The  corresponding  MTF  and  PSF  for  the 
VNIR  and  SWIR  multispectral  detectors  appear  in  Figures  8  through  11. 
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4.3  SCAN  MOTION  MTF 

Motion  of  the  satellite  as  it  collects  image  data  in  push-broom  mode  causes  a  slight  smearing  of 
the  image  in  the  in-track  direction.  The  geometric  footprint  of  each  pixel  on  the  ground  is 
convolved  with  a  vector  representing  the  motion  of  the  view  axis  during  the  integration  time  of  a 
single  frame.  The  MTF  of  this  motion  is  another  sine  function  in  the  in-track  direction  only: 

MTFij^t  -  sinc(av)  , 

where  v  is  the  angular  frequency  and  a  is  the  angle  subtended  by  the  ground  motion  vector. 

The  nominal  ground  speed  of  the  sub-satellite  point  is  6.75  km/s.  The  nominal  integration 
times  are  4  ms  and  1.33  ms  for  the  MS  and  Pan  bands,  respectively.  (Alternative  integration  times 
can  be  set  by  ground  command.)  The  nominal  motions  during  integration  are  thus  27.0  and  8.9 
meters,  respectively.  At  the  EO-1  altitude  of  705  km,  the  corresponding  values  of  a  are  38.28  and 
12.76  prad. 

4.4  UNCERTAINTIES  IN  THE  MODEL 

The  differences  between  the  modeled  MTF  and  the  measured  MTF  may  be  attributable  to  a 
number  of  causes.  In  this  section,  we  mention  several  areas  of  uncertainty  remaining  in  the  MTF 
model  after  we  have  adjusted  it  to  fit  the  measurements. 

4.4.1  Focus  Error 

The  necessity  of  setting  and  measuring  the  focus  of  the  instrument  in  a  clean  room  prevents  us 
from  testing  it  directly  by  observing  a  distant  object.  The  imaging  collimator  had  to  be  employed, 
and  the  estimated  error  in  setting  its  focus  is  0.145  pm  p-v,  over  the  12.5  cm  aperture. 

It  was  necessary  to  compensate  for  optical  power  induced  in  the  vacuum  chamber  window  when 
the  chamber  interior  was  cooled.’  Typically,  the  required  focus  correction  was  found  to  change  over 
the  course  of  five  hours  by  0.008  pm  p-v. 

The  focus  error  of  the  ALI  was  measured  for  one  location  on  the  FPA.  An  error  in  cutting  the 
focus  shim  at  the  correct  wedge  angle  (if  any)  would  cause  a  different  focus  error  at  other  points  on 
the  FPA. 

4.4.2  Collimator  Aberrations  and  Alignment 

We  lack  a  precise  description  of  the  random  aberrations  of  the  collimator  sufficient  to 
reconstruct  its  actual  WFE,  as  was  done  for  the  ALI.  Thus  those  aberrations  could  not  be 
incorporated  in  the  model  which  was  fitted  to  the  measured  MTFs.  The  spherical  aberration  of  the 
collimator  was  included,  but  it  increases  rapidly  as  the  collimator  exit  pupil  is  moved  away  from  its 
correct  registration  with  the  entrance  pupil  of  the  ALI.  At  the  time  most  of  these  measurements 
were  performed,  the  pupil  registration  was  uncertain  by  -^5  mm. 

4.4.3  Collimator  Ghosts 

Unwanted  reflections  at  surfaces  of  the  transmissive  elements  of  the  collimator  give  rise  to 
ghost  images.  The  most  noticeable  of  these  is  caused  by  light  reflected  from  the  back  face  of  the 
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beamsplitter,  instead  of  its  coated,  front  face.  This  displaced,  out-of-focus  ghost  has  an  intensity 
~1.2%  of  the  main  image.  We  believe,  without  doing  an  exhaustive  analysis,  that  these  ghost  images 
did  not  interfere  with  images  of  the  knife  edge  sufficiently  to  compromise  the  MTF  measurements. 

4.4.4  Scattered  Light  Effects 

Scattered  light  appears  to  have  reduced  the  main  portion  of  the  MTF  we  measure  by  as  much  as 
a  few  percent.  The  sources  of  this  scattering  are  the  ALI  mirrors,  as  mentioned  above,  and  a 
contamination  of  the  focal  plane.  The  exact  magnitude  of  the  scattering  fraction  as  well  as  the  true 
shape  of  the  scattering  MTF  at  angular  frequencies  less  than  ~  5  cycles/mrad,  is  somewhat  uncertain. 

4.4.4. 1  Mirrors 

The  knife-edge  scans  in  the  cross-track  direction  moved  the  image  of  the  edge  approximately 
1.2  mm  across  the  focal  plane  of  the  ALI.  This  is  not  sufficient  to  delineate  fully  the  effects  of  wide- 
angle  scattering  of  the  light.  Thus  the  reduction  of  the  specular  MTF  caused  by  light  scattering  at  the 
mirror  surfaces  was  underestimated  in  our  data  analysis. 

4. 4. 4.2  Contamination  on  filters 

Following  the  series  of  MTF  measurements  performed  during  December,  1998,  we  discovered 
that  an  unknown  contaminant  had  gradually  accumulated  on  the  top  surfaces  of  the  filters  which  are 
situated  just  above  the  FPA  detectors.  The  contamination  was  manifested  by  a  spatial  modulation  of 
the  light  from  the  internal  calibration  lamps,  amounting  to  several  percent.  We  obtained  images  of 
the  contaminant  by  placing  a  CCD  camera  at  the  focus  of  the  collimator  and  illuminating  the  FPA  of 
the  ALI  with  flashlights.  It  had  the  appearance  of  condensation  on  a  window  in  winter.  The 
modulation,  or  spatial  noise,  was  smaller  for  the  normal  //7.5  illumination  from  external  sources. 
Nevertheless,  it  is  probable  that  this  contaminant  caused  some  amount  of  wide-angle  scattering, 
which  because  of  its  proximity  to  the  detectors,  would  appear  to  extend  only  a  few  milliradians  in 
angle. 

The  contaminant  left  the  focal  plane  when  it  was  warmed  above  —5°  C.  Upon  re-cooling,  it 
reappeared,  but  not  so  strongly  as  before.  Evidence  of  similar  contamination  was  seen  during  thermal 
vacuum  tests  following  integration  of  the  ALI  with  the  spacecraft  bus.  It  is  not  known  where  the 
contaminant(s)  originates,  or  whether  it  will  reappear  on  orbit.  Provisions  have  been  made  to  warm 
the  focal  plane  periodically  on  orbit,  to  outgas  it  if  the  contamination  does  appear.  Thus  the  effects 
of  this  scattering  on  the  MTF  is  not  only  uncertain  by  a  few  percent,  but  could  be  variable  while  on 
orbit. 
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5  NUMERICAL  RESULTS 


5.1  WAVEFRONT  ERROR  OF  THE  ALI 

Procedures  written  in  IDL  are  used  to  construct  the  ALI  WFE  from  the  adjusted  set  of  Zernike 
coefficients  derived  from  interferograms  by  SSG.  From  the  WFE,  the  optical  point-spread  function 
(PSF)  and  MTF  of  the  ALI  telescope  are  obtained,  by  use  of  the  Fourier  transform  proeedure  in  IDL. 

A  surface  plot  of  the  WFE  near  the  middle  of  the  MS/Pan  array,  exeluding  tip,  tilt,  and  foeus,  is 
shown  in  Figure  12.  Contour  plots  of  the  WFE  at  the  Pan  arrays  of  SCAs  1  through  4  are  shown  for 
comparison  in  Figure  13.  There  are  small  differences  caused  by  the  fact  that  the  beam  footprint 
covers  different  parts  of  the  Ml  and  M3  mirrors  for  different  points  in  the  field  of  view. 


I 

i 
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Figure  12.  Surface  plot  of  the  wavefront  error  near  the  middle  of  the  MS/Pan  array. 
The  focus  term  is  not  included  here. 
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Figure  13.  Contour  plots  of  the  wavefront  error  at  the  centers  of  the  Pan  arrays.  Plots  a),  b),  c),  and  d)  represent 
SCAs  1,  2,  3,  and  4,  respectively.  The  interval  between  contours  is  63.3  nm  (1/10  HeNe  wave). 


When  the  modeled  MTF  was  compared  to  the  measured  Pan  MTFs  in  the  cross-track  and  in¬ 
track  directions,  an  apparent  inconsistency  was  foimd.  Much  better  fits  were  obtained  by  rotating  the 
WFE  array  by  90°.  On  consulting  SSG,  we  learned  that  the  orientation  of  the  interferogram  data  was 
ambiguous.  It  was  not  recorded,  because  the  interferograms  were  originally  taken  only  for  the  purpose 
of  obtaining  an  rms  wavefront  error  and  MTF  magnitude  at  a  specific  frequency.  Subsequently,  we 
used  the  rotated  WFE  for  our  MTF  model. 

The  optical  point-spread  function  of  the  Pan  band,  computed  from  the  WFE,  is  shown  in 
Figure  14.  The  real  part  of  the  corresponding  OTF  is  shown  in  Figure  15. 
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Figure  14.  Optical  point-spread  function  of  the  Pan  band,  computed  from  the  full  wavefront  error, 

including  the  focus  term. 


Figure  15.  Optical  transfer  function  of  the  Pan  band  (real  part),  computed  from  the  wavefront  error, 

including  the  focus  term. 


5.2  PANCHROMATIC  MTF 


Figure  1 6  shows  the  MTF  curves  derived  from  scans  of  the  four  Pan  arrays,  along  with  the 
average,  modeled  MTF  curve.  The  measured  MTF  curves,  and  corresponding  modeled  MTF  for  the 
four  individual  SCAs,  are  plotted  in  Figures  17  to  20. 
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Figure  16.  Measured  and  modeled  Average  MTFs  for  the  Pan  band,  and  measured  MTFs 

for  the  individual  SCAs. 
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Figure  1 7.  Measured  and  modeled  MTF  for  the  Pan  band,  SCA  1,  pixel  480. 
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Figure  18.  Measured  and  modeled  MTF  for  the  Pan  band,  SCA  2,  pixel  1440. 
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Figure  19.  Measured  and  modeled  MTF  for  the  Pan  band,  SCA  3,  pixel  2400. 
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Figure  20.  Measured  and  modeled  MTF  for  the  Pan  band,  SCA  4,  pixel  3360. 
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The  real  part  of  the  two-dimensional  Pan  MTF  at  the  center  of  the  array  computed  from  our 
model  is  shown  in  Figure  21.  The  small  imaginary  part  can  be  neglected,  given  the  uncertainties  in 
our  fitting  of  the  model  to  the  measurements.  The  real  part  of  the  MTF  is  symmetric  in  both  axes, 
the  imaginary  part  anti-symmetric.  MTF  values  over  one  quadrant  are  tabulated  at  reduced  resolution 
in  Appendix  A. 


Figure  21.  System  MTF  of  the  Panchromatic  band  near  the  center  of  the  detector  array. 


Figures  22  to  24  give  several  representations  of  the  Panchromatic  system  point-spread 
function  (PSF),  not  including  scan  motion  spreading. 
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Figure  22.  System  PSF  of  the  Panchromatic  band  near  the  center  of  the  detector  array. 
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Figure  23.  System  PSF  of  the  Panchromatic  band  near  the  center  of  the  detector  array. 
The  contour  levels  are  at  5,  10,  30,  50,  70,  90  and  95%  of  the  peak. 


Figure  24.  System  PSF  of  the  Panchromatic  band  near  the  center  of  the  detector  array. 


5.3  MULHSPECTRALMTFs 


Summary  plots,  similar  to  Figure  16,  of  the  MTFs  of  each  of  the  multispectral  bands  are  shown 
in  Figures  25  to  33.  The  real  part  of  the  two-dimensional  average  MTFs  for  bands  1’,  4’,  and  7  are 
plotted  in  Figures  34  to  36.  Appendix  A  contains  the  tabulated  values. 


Cross— track  MTF 


In  —  trock  MTF 


Figure  25.  Measured  and  modeled  Average  MTFs  for  band  1  and  measured  MT’s  for  the  individual  SCAs. 
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Figure  26.  Measured  and  modeled  Average  MTFs  for  band  1,  and  measured  MTFs  for  the  individual  SCAs. 
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Figure  27.  Measured  and  modeled  Average  MTFs  for  band  2,  and  measured  MTFs  for  the  individual  SCAs. 
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Figure  29.  Measured  and  modeled  Average  MTFs  for  band  4,  and  measured  MTFs  for  the  individual  SCAs. 
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Figure  32.  Measured  and  modeled  Average  MTFs  for  band  5,  and  measured  MTFs  for  the  individual  SCAs. 
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Figure  33.  Measured  and  modeled  Average  MTFs  for  band  7,  and  measured  MTFs  for  the  individual  SCAs. 
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Figure  38.  System  PSF  of  band  I  ’  near  the  center  of  the  detector  array. 
The  contour  levels  are  at  5,  10,  30,  50,  70,  90  and  95%  of  the  peak. 


Cross— Track  (/Am) 

Figure  39.  System  PSF  of  band  1  ’  near  the  center  of  the  detector  array.  This  cross-section  is  through  center 

in  the  cross-track  direction. 
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Figure  42.  System  PSF  of  band  4  ’  near  the  center  of  the  detector  array. 
This  cross-section  is  through  center  in  the  cross-track  direction. 


Figure  43.  System  PSF  of  band  7  (2.2  fim)  near  the  center  of  the  detector  array. 
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Figure  44.  System  PSF  of  band  7  near  the  center  of  the  detector  array. 
The  contour  levels  are  at  5,  10,  30,  50,  70,  90  and  95%  of  the  peak. 
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Figure  45.  System  PSF  of  band  7  near  the  center  of  the  detector  array. 
This  cross-section  is  through  center  in  the  cross-track  direction. 


6  REFERENCES 


1.  Bicknell,  W.  E.,  Digenis,  C.  J.,  Forman,  S.  E.,  and  Lencioni,  D.  E.,  EO-1  Advanced  Land  Imager, 

Proc.  SPIE,  vol.  3750,  pp.  80-88,  July  1999. 

2.  Lencioni,  D.  E.,  Hearn,  D.  R.,  Mendenhall,  J.  A.,  and  Bicknell,  W.  E.,  EO-1  Advanced  Land 

Imager  Calibration  and  Performance  Overview,  Proc.  SPIE,  vol.  3750,  pp.  89-96,  July  1999. 

3.  Hearn,  D.  R.,  Mendenhall,  J.  A.,  and  Willard,  B.  C.,  Spatial  Calibration  of  the  EO-1  Advanced 

Land  Imager,  Proc.  SPIE,  vol.  3750,  pp.  97-108,  July  1999. 

4.  Willard,  B.  C.,  Wide  field-of-view  Schmidt-sphere  imaging  collimator,  Proc.  SPIE,  vol.  3750,  pp 

286-296,  July  1999. 

5.  Lambda  Research  Corporation,  EO-1  Stray  Light  Analysis  Report  No.  2,  Feb.  20,  1998. 

6.  Tamowski,  G.  J.,  and  Maschhoff,  K.  R.,  FPA/Dewar  MTF  Component  Analysis,  Proc.  IRIS 

Detectors,  pp.  69-94,  1996. 

7.  Heam,  D.  R.,  Vacuum  window  optical  power  induced  by  temperature  gradients,  Proc.  SPIE,  vol. 

3750,  pp  297-308,  July  1999. 


39 


APPENDIX  A 


SPATIAL  TRANSFER  FUNCTIONS  OF  THE 
EO-1  ADVANCED  LAND  IMAGER 

The  spatial  transfer  functions  (STFs)  of  the  ALI  are  tabulated  on  the  following  pages. 
These  nmnerical  values  are  obtained  with  our  analytical  model,  as  described  in  Section  4  of 
this  report.  The  location  on  the  focal  plane  that  was  used  for  these  calculations  is  SCA  2, 
pixel  635,  which  is  in  the  center  of  the  cross-track  extent  of  the  populated  part  of  the  focal 
plane.  The  STF  varies  only  slightly  over  the  field,  so  the  values  given  here  can  be  considered 
fully  representative. 

Tables  A-1  to  A-10  list  the  STFs  of  each  of  the  bands.  The  complex  STF  is  given  for 
one  quadrant  of  frequency  space.  The  real  part  is  symmetric,  the  imaginary  part  anti¬ 
symmetric  in  both  axes.  The  frequencies  listed  in  the  top  rows  and  left  columns  of  the  tables 
are  first,  spatial  frequency  at  the  focal  plane  (in  cycles/millimeter),  and  second,  angular 
frequency  in  object  space  (in  cycles/milliradian).  (They  are  related  by  the  focal  length  of  the 
instrument,  0.946  meters.)  The  frequency  intervals  are  one-fourth  of  the  Nyquist  frequency 
in  each  case. 

In  the  main  body  of  each  table,  the  STF  value  is  listed  with  the  real  part  above,  and  the 
imaginary  part  below  in  the  same  cell. 
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APPENDIX  B 


ESTIMATED  ZERNIKE  POLYNOMIAL  COEFTTCIENTS  OF  WAVEFRONT  ERRORS 

SSG,  Inc.  used  a  laser  unequal-path  interferometer  (LUPI)  to  characterize  the  WFE  of  the  ALI 
after  final  alignment  of  the  optical  system.  They  obtained  interferograms  at  eleven  points  distributed 
over  the  1.26°xl5°  field  of  view  of  the  optical  system.  The  interferograms  were  referenced  to  a 
fixture  representing  the  focal  plane.  The  WFE  derived  from  each  interferogram  was  fitted  with  a  set 
of  37  Zemike  polynomials. 

The  coordinates  on  the  focal  plane  where  the  interferograms  were  made  are  given  in  Table  B-1. 
Those  positions  are  refeired  to  the  center  of  the  15°xl.26°  field,  with  Xsc  tke  in-track  direction,  and 
the  SCAs  at  the  +Ysc  end  of  the  focal  plane.  Table  B-2  lists  the  Zemike  polynomial  definitions  used 
here.  The  derived  coefficients  for  Zemike  terms  4  through  37  are  shown  in  Table  B-3.  (Terms  1,  2, 
and  3  are  set  to  zero,  since  they  only  represent  average  phase  and  tilt  of  the  wavefront.) 


TABLE  B-1 

Measurement  Positions  in  the  Focal  Plane 


Position 

Xsc  (nim) 

Ysc  (mm) 

FP  1 

-10.16 

127.00 

FP2 

0.00 

127.00 

FP3 

10.16 

127.00 

FP4 

-10.16 

63.50 

FP5 

0.00 

63.50 

FP6 

10.16 

63.50 

FP7 

-10.16 

0.00 

FP8 

0.00 

0.00 

FP9 

10.16 

0.00 

FP  10 

-10.16 

-127.00 

FP  12 

10.16 

-127.00 
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TABLE  B-2 

Zernike  Polynomial  Definitions 


Term 

Polynomial 

Term 

Polynomial 

1 

2 

1 

rcos(0) 

14 

1 

lOr^-12r^+3r]cos(0) 

lOr^-12r^+3r)sin(0) 

20r^-30/+12r2-l 

3 

4 

rsin(0) 

2r^-l 

■ 

cos(20) 

17 

cos(40) 

6 

r“  sin(20) 

18 

sin(40) 

7 

^3r^  -2rjcos(0) 

19 

V-4r2)cos(30) 

8 

^3r^  -2rjsin(0) 

20 

V-4r^)sin(30) 

9 

6r^-6r^+l 

21 

15/-2Or'^  +  6r2)cos(20) 

10 

cos(30) 

22 

'l5r^-20r'^+6r~)sin(2e) 

11 

sin(30) 

23 

^35r^  -  60r^  +  30r^  -  4rjcos(0) 

12 

-3r^jcos(20) 

24 

i  35r’  -  60r^  +  30r^  -  4r)sin(0) 

13 

(4r^-3r2)sin(20) 

25 

70r^  -  140r^  +  90/  -  20r^  + 1 

cos(50) 

sin(50) 

(6r^-5r^)cos(40) 
l6r^  -  5r^jsin(40) 

(21r^-30r^  +  10r^)cos(36>) 

(21r^-3Or5  +  lOr^)sin(30) 

(56r^  -  105r^  +  60r^  -  10r^)cos(26>) 

(56r^  -  105r^  +  60r'^  -  10r^)sin(26l) 

(l26r^  -  280r’  +  210r^  -  60^^  +  5r)cos(0) 

(l26r^  -  280r^  +  210r^  -  60r^  +  5r)sin(0) 

252r'®  -  630r^  +  560r^  -  210/  +  30r^  - 1 

924r^^  -  2772r^®  +  3150r^  -  1680r^  +  420r'^  -  42r^  + 1 
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TABLE  B-3 


Estimated  Zernike  Polynomial  Coefficients 


Term 

FP1 

FP2 

Pos 

FP3 

ition 

FP4 

FP5 

FP6 

4 

-0.7885 

-0.7765 

-0.7463 

-0.7080 

-0.7567 

-0.7572 

5 

0.1889 

0.1747 

0.1559 

-0.0811 

-0.0642 

-0.0721 

6 

-0.0999 

-0.0197 

0.0176 

0.0962 

0.1117 

0.0993 

7 

0.1275 

0.1197 

0.1179 

0.0164 

-0.0023 

-0.0269 

8 

0.0573 

0.0019 

-0.0002 

0.0284 

0.0150 

-0.0228 

9 

0.0639 

0.0769 

0.0746 

-0.0022 

0.0126 

0.0197 

10 

0.1184 

0.1021 

0.0615 

-0.0934 

-0.1030 

-0.0813 

1 1 

-0.0217 

0.0163 

0.0093 

0.0373 

0.0470 

0.0696 

12 

0.1367 

0.1438 

0.1448 

0.0472 

0.0201 

-0.0010 

13 

-0.0247 

0.0110 

0.0593 

0.0188 

0.0466 

0.0701 

14 

0.0468 

0.0174 

0.0176 

-0.0322 

-0.0406 

-0.0376 

15 

-0.0093 

-0.0454 

0.0176 

0.0189 

0.0218 

0.0175 

16 

0.0373 

0.0023 

-0.0049 

0.0287 

0.0282 

0.0335 

17 

0.1495 

0.1184 

0.0712 

-0.0075 

0.0520 

0.0986 

18 

0.0055 

0.0098 

-0.0339 

0.0323 

0.0095 

0.0049 

19 

-0.0355 

0.0649 

0.0797 

0.0157 

0.0205 

0.0122 

20 

-0.0230 

0.0202 

0.0147 

0.0098 

0.0269 

0.0195 

21 

-0.0355 

0.0190 

0.0655 

0.0054 

-0.0009 

-0.0054 

22 

0.0413 

0.0097 

0.0018 

0.0011 

-0.0116 

-0.0089 

23 

-0.0160 

-0.0113 

0.0133 

0.0217 

0.0102 

0.0102 

24 

0.0452 

0.0400 

0.0182 

-0.0097 

-0.0182 

-0.0199 

25 

-0.0018 

0.0012 

0.0073 

0.0047 

-0.0038 

-0.0144 

26 

0.0372 

-0.0232 

-0.0386 

-0.0133 

-0.0235 

0.0031 

27 

0.0931 

0.0793 

0.0476 

-0.0287 

-0.0592 

-0.0445 

28 

0.0140 

0.0506 

0.0227 

0.0050 

-0.0109 

-0.0421 

29 

-0.0397 

-0.0083 

0.0181 

0.0178 

0.0104 

-0.0058 

30 

-0.0554 

-0.0006 

0.0211 

-0.0129 

-0.0024 

-0.0114 

31 

0.0200 

-0.0097 

-0.0224 

-0.0168 

-0.0091 

-0.0162 

32 

-0.0367 

-0.0014 

0.0226 

-0.0167 

0.0001 

0.0076 

33 

0.0332 

0.0292 

-0.0032 

0.0288 

0.0129 

0.0044 

34 

-0.0710 

-0.0423 

-0.0120 

-0.0204 

-0.0049 

0.0023 

35 

-0.0494 

-0.0053 

-0.0262 

-0.0179 

-0.0078 

-0.0166 

36 

-0.0103 

-0.0279 

0.0007 

-0.0076 

-0.0066 

0.0043 

37 

0.0074 

0.0195 

0.0058 

-0.0001 

0.0006 

-0.0139 

(continued) 
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TABLE  B-3 

Estimated  Zernike  Polynomial  Coefficients  (continued) 


Position 

Term 

FP7 

FP8 

1  FP9 

FP  10 

FP  12 

-0.5369 

-0.0484 

0.0171 

-0.0316 

-0.0264 

-0.0489 

-0.0442 

-0.0752 

0.0198 

0.0712 

0.0290 

-0.0122 

0.0643 

-0.0522 

0.0376 

0.0064 

0.0239 

0.0042 

-0.0359 

-0.0083 

-0.0108 

-0.0004 

0.0289 

0.0117 

-0.0165 

0.0059 

-0.0103 

-0.0094 

-0.0239 

0.0079 

-0.0094 

-0.0299 

-0.0045 

-0.0013 


-0.5140 

-0.0134 

0.0186 

-0.0535 

-0.0469 

-0.0443 

-0.0369 

-0.0906 

-0.0164 

0.0717 

0.0196 

-0.0303 

0.0630 

-0.0156 

0.0209 

0.0017 

0.0394 

0.0100 

-0.0171 

-0.0028 

-0.0062 

-0.0131 

-0.0018 

0.0029 

-0.0078 

-0.0179 

-0.0081 

0.0004 

-0.0039 

-0.0104 

-0.0216 

-0.0141 

-0.0173 

0.0067 


-0.5408 

0.0090 

0.0180 

-0.0554 

-0.0925 

-0.0328 

-0.0184 

-0.0851 

-0.0507 

0.0587 

0.0029 

-0.0247 

0.0579 

0.0230 

0.0107 

0.0001 

0.0433 

0.0153 

-0.0177 

0.0053 

0.0192 

-0.0125 

-0.0249 

0.0021 

-0.0006 

-0.0240 

-0.0131 

-0.0071 

0.0094 

-0.0158 

-0.0174 

-0.0082 

-0.0114 

0.0065 


-0.3554 

0.1974 

-0.0154 

-0.1158 

0.1282 

0.0463 

-0.0219 

0.0581 

0.0773 

0.0112 

0.0400 

0.0097 

0.0332 

0.1109 

-0.0263 

0.0287 

-0.0266 

-0.0075 

-0.0352 

-0.0102 

-0.0142 

-0.0155 

-0.0076 

0.1166 

-0.0726 

0.0018 

-0.0208 

-0.0080 

-0.0012 

0.0531 

0.0026 

-0.0514 

0.0048 

0.0268 


-0.3121 

0.2517 

-0.0718 

-0.1316 

0.0171 

0.0460 

0.0051 

0.0623 

0.0699 

-0.0616 

0.0525 

-0.0134 

0.0272 

0.0636 

-0.0509 

0.0175 

0.0215 

0.0171 

-0.0119 

-0.0205 

0.0208 

0.0017 

0.0485 

0.0532 

0.0184 

0.0010 

-0.0107 

-0.0050 

-0.0264 

-0.0296 

-0.0003 

0.0006 

0.0029 

-0.0044 
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